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We study the shear and bulk viscosities of partonic and hadronic matter as functions of tempera- 
ture T within the Parton-Hadron-String Dynamics (PHSD) off-shell transport approach. Dynamical 
hadronic and partonic systems in equilibrium are studied by the PHSD simulations in a finite box 
with periodic boundary conditions. The ratio of the shear viscosity to entropy density n(T) / s(T) 
from PHSD shows a minimum (with a value of about 0.1) close to the critical temperature T c , while 
it approaches the perturbative QCD (pQCD) limit at higher temperatures in line with lattice QCD 
results. For T < T c , i.e. in the hadronic phase, the ratio rj/s rises fast with decreasing temperature 
due to a lower interaction rate of the hadronic system and a significantly smaller number of degrees- 
of-freedom. Within statistics, we obtain practically the same results in the Kubo formalism and 
in the relaxation time approximation. The bulk viscosity £(T) - evaluated in the relaxation time 
approach - is found to strongly depend on the effects of mean fields (or potentials) in the partonic 
phase. We find a significant rise of the ratio £(T)/s(T) in the vicinity of the critical temperature 
T c , when consistently including the scalar mean-field from PHSD, which is also in agreement with 
that from 1QCD calculations. Furthermore, we present the results for the ratio (n + 3(/4)/s, which 
is found to depend non-trivially on temperature and to generally agree with the 1QCD calculations 
as well. Within the PHSD calculations, the strong maximum of (,{T)/n(T) close to T c has to be 
attributed to mean-fields (or potential) effects that in PHSD are encoded in the temperature de- 
pendence of the quasiparticle masses, which is related to the infrared enhancement of the resummcd 
(effective) coupling g(T). 



I. INTRODUCTION 



High energy heavy-ion reactions are studied experi- 
mentally and theoretically to obtain information about 
the properties of nuclear matter under the extreme condi- 
tions of high baryon density and/or temperature. Ultra- 
rclativistic heavy-ion collisions at the Relativistic Heavy- 
Ion Collider (RHIC) and the Large Hadron Collider 
(LHC) at CERN have produced a new state of matter, 
the quark-gluon plasma (QGP) , for a couple of fm/c. The 
produced QGP shows features of a strongly-interacting 
fluid unlike a weakly-interacting parton gas [H-Q • Large 
values of the observed azimuthal asymmetry of charged 
particles in momentum space, i.e. the elliptic flow V2 
[5j-l9|, could quantitatively be well described by hydro- 



dynamics up to transverse momenta on the order of 1.5 
GeV/c [10l4l5j . A perfect fluid has been defined as hav- 
ing a zero shear viscosity 77; yet semiclassical arguments 
have been given suggesting that the shear viscosity can- 
not be arbitrarily small [16]. Indeed, the lower bound for 
the shear viscosity to entropy density ratio rj/s = 1 /An 
was conjectured by Kovtun-Son-Starinets (KSS) [l?} for 
infinitely coupled supersymmetric Yang-Mills gauge the- 
ory based on the AdS/CFT duality conjecture. On the 
basis of holographically dual computations [l8[ , also for 
the bulk viscosity of strongly-coupled gauge theory plas- 
mas a lower bound was conjectured: (/rj > 2(1/3 — c? s ), 
where c s is the speed of sound. Empirically, relativis- 
tic viscous hydrodynamic calculations - using the Israel- 
Stewart framework - require a very small but finite rj/s 
of 0.08 — 0.24 in order to reproduce the RHIC elliptic 
flow V2 data [l9l - [22| . The main uncertainty in these es- 
timates results from the equation of state and the initial 
conditions employed in the hydrodynamical calculations 
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as well as in the temperature dependence of rj/s. 

Thus not only the absolute values of the shear and bulk 
viscosities are of great interest, but also the temperature 
dependence of these coefficients, which is expected to be 
quite strong. There is evidence from atomic and molec- 
ular systems that rj/s should have a minimum in the 
vicinity of the phase transition or - in case of strong in- 
teractions at vanishing chemical potential - of the rapid 
crossover between hadronic matter and the quark-gluon 
plasma [23T - [25l |. Furthermore, it is argued that the ratio 
of the bulk viscosity to entropy density £/s should have 
a maximum close to T c - as suggested by lattice QCD - 
and might even diverge in case of a second-order phase 
transition 26-31J. Such a peak in the bulk viscosity can 
lead to instabilities in viscous hydrodynamics simulations 
for heavy-ion collisions and possibly to clusterization ef- 
fects [If. 

Shear and bulk viscosities of strongly-interacting sys- 
tems have been evaluated within different approaches. 
Calculations have been performed at high temperatures, 
where perturbation theory can be applied [33l |34| . as 
well as at extremely low temperatures [3J-|36|. First re- 
sults for shear and bulk viscosities obtained within lat- 
tice QCD (1QCD) simulations above the critical tem- 
perature of pure gluon matter have been presented in 
[37H40l |. There are several methods for the calculation 
of shear and bulk viscosities for strongly-interacting sys- 
tems: the Relaxation Time Approximation (RTA) [42j . 
the Chapmann-Enskog (CE) method [43| and the Green- 
Kubo approach [4J, |45[ . The RTA method has been used 
to calculate the viscosity fla . [30 l KMa - l5l| as well as the 
Green-Kubo approach [J, [421 [52- 56] for both hadronic 
and partonic matter providing a rough picture of the 
transport properties of strongly-interacting matter. 

In this study we will calculate the shear and bulk vis- 
cosities as a function of temperature (or energy density) 
with the Parton-Hadron- String Dynamics (PHSD) trans- 
port approach that has provided a good description of 
collective flow properties and differential particle spec- 
tra in nucleus- nucleus collisions from lower SPS to RHIC 
energies [57l - [6l| . In this approach the shear and bulk 
viscosities do not enter as external parameters but are 
generic properties of the matter under consideration and 
can be calculated for systems in equilibrium as a func- 
tion of temperature explicitly without incorporating any 
additional parameters. 

The paper is organized as follows. In Sec. II we provide 
a brief reminder of the off-shell dynamics and the ingre- 
dients of the PHSD transport approach. We then first 
present in Sec. Ill the actual results for the shear and 
bulk viscosities in 'infinite' parton-hadron matter within 
the PHSD employing the Green-Kubo formalism and the 
Relaxation Time Approximation (RTA) and compare to 
available lattice QCD results. The summary and conclu- 
sions are given in Sec. IV. 



II. THE PARTON-HADRON-STRING 
DYNAMICS (PHSD) TRANSPORT APPROACH 

In this work we extract the shear and bulk viscosities 
for 'infinite' parton-hadron matter employing different 
methods within the PHSD transport approach [5?], HU , 
which is based on generalized transport equations on the 
basis of the off-shell Kadanoff-Baym equations [HI, [63| 
for Green's functions in phase-space representation (in 
the first order gradient expansion, beyond the quasi- 
particle approximation). The approach consistently de- 
scribes the full evolution of a relativistic heavy-ion colli- 
sion from the initial hard scatterings and string formation 
through the dynamical deconfinement phase transition 
to the strongly-interacting quark-gluon plasma (sQGP) 
as well as hadronization and the subsequent interactions 
in the expanding hadronic phase. In the hadronic sec- 
tor PHSD is equivalent to the Hadron-String-Dynamics 
(HSD) transport approach !64s,.6H _ a covariant extension 
of the Boltzmann-Uehling-Uhlenbeck (BUU) approach 
[o6| - that has been used for the description of pA and 
AA collisions from SIS to RHIC energies in the past. In 
PHSD the partonic dynamics is based on the Dynamical 
QuasiParticle Model (DQPM) [67M69j. which describes 
QCD properties in terms of single-particle Green's func- 
tions (in the sense of a two-particle irreducible (2 PI) 
approach) and reproduces lattice QCD results - includ- 
ing the partonic equation of state - in thermodynamic 
equilibrium. 



A. Reminder of the DQPM 

In the scope of the DQPM the running coupling con- 
stant g 2 (squared) for partons is approximated (for T > 
Tc) by 
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where the parameters A = 2.42 and T s /T c — 0.56 have 
been extracted from a fit to the lattice data for purely 
gluonic systems (Nf=0) as described in Ref. [69j. In 
Eq. (fT]), N c = 3 stands for the number of colors, T c is 
the critical temperature (= 158 MeV), while Nf denotes 
the number of flavors. In the actual PHSD calculations 
for Nf = 3 we employ a slightly different analytical form 
for g 2 (T/T c ) that has been fitted to the lattice data from 
Ref. ■ For the details we refer the reader to Ref. [7(j . 

The functional forms for the dynamical quasiparticle 
masses (for gluons and quarks) are chosen so that they 
become identical to the perturbative thermal masses in 
the asymptotic high-temperature regime, i.e., for gluons 
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and for quarks (antiquarks) 
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but the running coupling g(T/T c ) is the resummed cou- 
pling of Eq. (1). The effective quarks, antiquarks and 
gluons in the DQPM have finite widths, which for p q = 
are adopted in the form [71 1 



r 9 (T) 



1 Ar g a T (2c 
-N P ^— In -s 



3 8tt \g 
IN 2 -I g 2 T 
3 27V C 8T 



1 



(4) 
(5) 



where the parameter c = 14.4 is related to a magnetic 
cut-off (see Ref. 0]). 

In line with Ref. [69| , the parton spectral functions are 
no longer 5-functions in the invariant mass squared but 
have a Breit-Wigner form 
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index j stands for quarks, antiquarks and gluons (j — 
q,q,g). The spectral function ([5]) is antisymmetric in uj 
and normalized as 
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The parameters Tj , My from the DQPM have been de- 
fined above. Note, however, that the decomposition of 
the total width Tj into the collisional width (due to elas- 
tic and inelastic collisions) and the decay width is not 
addressed in the DQPM. The effective cross sections for 
each of the various partonic channels as a function of the 
energy density e, which fixes the partial widths of the dy- 
namical quasiparticles as well as the various interaction 
rates, has been determined in Ref. 741. 



B. Hadronization in PHSD 

The hadronization, i.e. the transition from partonic to 
hadronic degrees-of-freedom and vice versa is described 
in PHSD by covariant transition rates for the fusion 
of quark- antiquark pairs to mesonic resonances or three 
quarks (antiquarks) to baryonic states [58j . e.g. for q + q 
fusion to a meson m of four-momentum p = (w,p) at 



space-time point x = (t, x 
dN m (x,p) 
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xN q (xq,Pq)Nq(xq,Pq)S(&avor, color) . (8) 
In ([3} we have introduced the shorthand notation, 
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where ■ denotes a summation over discrete quantum 
numbers (spin, flavor, color); Nj(x,p) is the phase-space 
density of parton j at space-time position x and four- 
momentum p. In © £ (flavor, color) stands symbolically 
for the conservation of flavor quantum numbers as well 
as color neutrality of the formed hadron m, which can be 
viewed as a color-dipole or 'pre-hadron'. Furthermore, 
Vq q (p p ) is the effective quark- antiquark interaction from 
the DQPM (displayed in Fig. 10 of Ref. [68)]) as a func- 
tion of the local parton (q + q + g) density p p (or en- 
ergy density). Furthermore, W m (x,p) is the dimension- 
less phase-space distribution of the formed 'pre-hadron', 
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with £ = xi - x 2 = x q - x q and p$ = (pi - Pp)/2 = 
(p q —pq)/2 (which had been introduced in Ref. [72j). The 
width parameter b has been fixed by \J (r 2 ) = b = 0.66 fm 
(in the rest frame) , which corresponds to an average rms 
radius of mesons. We note that the expression (fT0|) cor- 
responds to the limit of independent harmonic oscillator 
states and that the final hadron-formation rates are ap- 
proximately independent of the parameter b within rea- 
sonable variations. By construction the quantity (|10p is 
Lorentz invariant; in the limit of instantaneous 'hadron 
formation', i.e. £° = 0, it provides a Gaussian drop- 
ping in the relative distance squared (ri — r2) 2 . The 
four-momentum dependence reads explicitly (except for 
a factor 1/2) 

{E l - E 2 f - (pi - p 2 ) 2 - (Mi - M 2 ) 2 < (11) 

and leads to a negative argument of the second expo- 
nential in (|10p favoring the fusion of partons with low 
relative momenta p q — pq = p\ — p 2 . 

Note that due to the off-shell nature of both partons 
and hadrons, the hadronization process obeys all conser- 
vation laws (i.e. the four-momentum conservation and 
the flavor current conservation) in each event, the de- 
tailed balance relations and the increase in the total en- 
tropy S for rapidly expanding systems. The physics be- 
hind (O is that the inverse reaction, i.e. the dissolution 
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FIG. 1: (color online) The correlation functions 
(-K xy (0)-K xv {t)), which are normalized by (ir xy (0) 2 ), as 
a function of time from the PHSD simulations in the box 
(open symbols) for systems at different energy densities. The 
corresponding exponential fits are given by dashed lines; the 
extracted relaxation times r are given, too. 



of hadronic states to quark-antiquark pairs (in case of 
mesons), at low energy density is inhibited by the large 
masses of the partonic quasiparticles according to the 
DQPM. Vice versa the resonant q — q pairs have a large 
phase-space to decay to several 0~ octet mesons. Wc 
recall that the transition matrix element becomes huge 
below the critical energy density [58]. For further details 
on the PHSD off-shell transport approach and hadroniza- 
tion we refer the reader to Refs. |H3, HI, [70l.l73l.l74l. 



III. CALCULATION OF SHEAR AND BULK 
VISCOSITY COEFFICIENTS 
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FIG. 2: (color online) The relaxation time r as a function of 
the number of test-particles obtained by the PHSD simula- 
tions in the box (symbols) for systems at different energies 
density. The dashed lines provide the convergent values for 
the relaxation time r. 



where the particle abundances, spectra, fluctuations and 
spectral functions were studied. In the present work, we 
extend our investigations to the calculation of transport 
coefficients. 



A. The Kubo formalism 

The Kubo formalism relates linear transport coeffi- 
cients such as heat conductivity, shear and bulk viscosi- 
ties to non-equilibrium correlations of the corresponding 
dissipative fluxes and treats dissipative fluxes as pertur- 
bations to local thermal equilibrium pi . [45| . The Green- 
Kubo formula for the shear viscosity r] is as follows [75| : 



In this Section we concentrate on the extraction of 
the shear and bulk viscosities for 'infinite' parton-hadron 
matter employing the Green-Kubo formalism and the re- 
laxation time approximation. We simulate the 'infinite' 
matter within a cubic box with periodic boundary con- 
ditions at various values for the energy density within 
PHSD. The size of the box is fixed to 9 3 fm 3 . The initial- 
ization is done by populating the box with light (u, d) and 
strange (s) quarks, antiquarks and gluons. If the energy 
density in the system is below the critical energy density 
(e c sa 0.5 GeV/fm 3 ), the evolution proceeds through the 
dynamical phase transition (as described in Section III B[) 
and ends up in an ensemble of hadrons. The system is 
initialized slightly out of equilibrium and, at all energy- 
densities, approaches kinetic and chemical equilibrium 
during it's evolution within PHSD as was shown in our 
previous investigations in Ref. [741 ] . After equilibration, 
the properties of the system at given temperature can be 
studied. For more details we refer the reader to Ref. [7J| , 
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where T is the temperature of the system and (...) de- 
notes the ensemble average in thermal equilibrium. In 
Eq. p^|) Tr xy is the shear component (non-diagonal spa- 
cial part) of the energy momentum tensor n^: 
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where the scalar mean-field U s (from PHSD) enters in 
the energy E = ^/p 2 + U%. 

In our numerical simulation - within test-particles rep- 
resentation - the volume averaged shear component of 
the energy momentum tensor can be written as 
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FIG. 3: (color online) The shear viscosity rj as a function of 
temperature from the PHSD simulations in the box for vari- 
ous volumes of the box: V = 125 GeV/fm 3 (green squares), 
V = 343 GeV/fm 3 (red circles) and V = 729 GeV/fm 3 (blue 
triangles) . 



{Ntest = 500), which leads to reliable (within statistical 
error bars) results. 

We also note that our numerical results for rj do not 
depend on the volume V of the box within reasonable 
variations by factors of 6 as shown in Fig. [3] 



B. The relaxation time approximation 

The starting hypothesis of the relaxation time approx- 
imation (RTA) is that the collision integral can be ap- 
proximated by 



cm = - 
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where r is the relaxation time. In this approach it has 
been shown that the shear and bulk viscosities (without 
mean-field or potential effects) can be written as [76|-[zH : 
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where V is the volume of the system and the sum is over 
all particles in the box at time t. Note that the scalar 
mean-field contribution U s only enters via the energy E. 
Taking into account that point particles are uniformly 
distributed in our box (implying ir xy (r, t) — n xy (t))), we 
can simplify the Kubo formula for the shear viscosity to 

oo 

V = ¥.Jdt{n**(0)**y(t)). (15) 
o 

The correlation functions (ir xy (0)ir xy (t)) are empirically 
found to decay almost exponentially in time 

(7^(0)7^)) = (n xy (0)ir xy (0)) e~ t/T (16) 

as shown in Fig. [TJ where r is the respective relaxation 
times. Finally, we end up with the Green-Kubo formula 
for the shear viscosity: 

V = ^ xy (0f)r, (17) 

which we use to extract the shear viscosity from the 
PHSD simulations in the box at given energy density. 
Note that the temperature T is uniquely related to the 
energy density e(T) in PHSD (in thermodynamic equi- 
librium) . 

We check the numeric stability of the method by plot- 
ting the respective relaxation times t, extracted from 
the PHSD simulations in the box, as a function of the 
number of test-particlcs in Fig. [2] The results for the re- 
laxation time r converge for N test ^ 400 independent of 
the energy density. In this study, we use a high amount 
of microcanonical simulations in the ensemble average 



where the sum is over particles of different type a (in our 
case, a = q,q,g). In the PHSD transport approach the 
relaxation time is given by: 

r a (T)=T-\T), (21) 

where T a (T) is the width of particles of type a = q,q,g 
as defined by (jHEJ)- In °u r numerical simulation - within 
test-particles representation - the volume averaged shear 
and bulk viscosities are given by the following expres- 
sions: 

1 W i 14 
l—l 1 

i N r- 1 

c-r^E^rtd- 3 ^ 2 -^] 2 ' ( 23 ) 

i— 1 % 

where the speed of sound v s = v s (T) is taken from lat- 
tice QCD [zU or the DQPM alternatively. Note that 
v s (T) from both approaches is practically identical since 
it is governed by the DQPM, which reproduces the lattice 
QCD results. 

C. Results for the shear viscosity 

In Fig. U we present the shear viscosity to entropy 
density ratio rj/s as a function of temperature T of 
the system extracted from the PHSD simulations in 
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FIG. 4: (color online) The shear viscosity to entropy den- 
sity ratio t)/s as a function of temperature of the system ob- 
tained by the PHSD simulations using different methods: the 
relaxation time approximation (red line+diamonds) and the 
Kubo formalism (blue line+dots). The others symbols de- 
note lattice QCD data for pure SU C (3) gauge theory from 
[37| (magenta squares), from [3!| (open and full triangles), 
and from [40l ] (black open and full circles). The orange 
dashed line demonstrates the Kovtun-Son-Starinets bound 
[13 {v/ s )kss = l/(47r). For comparison, the results from 
the virial expansion approach (green line) [2J] are shown as a 
function of temperature, too. 



the box employing the relaxation time approximation 
(red line+diamonds) and the Kubo formalism (blue 
line+dots). For comparison, the results from the virial 
expansion approach [24| are given by the green line as 
well as lattice QCD data for pure SU C (3) gauge theory. 
We find that all approaches give roughly the same r]/s as 
a function of temperature within error bars. The results 
for T < T c stem from PHSD in the Kubo framework 
and rapidly rise with decreasing temperature since the 
hadronic system is less strongly-interacting and the en- 
tropy density is lower, too, due to the lower number of 
degrees-of-freedom in the hadronic phase. 



FIG. 5: (color online) The bulk viscosity to entropy den- 
sity ratio as a function of temperature T of the system ex- 
tracted from the PHSD simulations in the box using the re- 
laxation time approximation including mean-field effects (red 
line+diamonds) and without potential effects (blue line+open 
triangles) . The available lattice QCD data from [3|| are given 
by green squares and from [401 ] by black open and full circles, 
respectively. 



with potential effects reads 
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In the numerical simulation the volume averaged bulk 
viscosity (including the mean-field effects from PHSD) is 
evaluated as 
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By using the DQPM expressions for the masses of quarks 
and gluons (for [i q = 0) 

(T/T c ) = \g 2 {T/T c )T 2 , m 2 (T/T c ) = -a 2 {T/T c )T 2 



D. Mean-field or potential effects 



We recall that for vanishing quark chemical potential 
the partonic mean-fields are essentially of scalar type and 
vector or tensor fields are subleading since the average 
quark current is zero. Furthermore, partonic mean-fields 
affect the bulk viscosity but not the shear viscosity (ex- 
cept for a contribution in the energy E in the dominator). 
According to [78|, the expression for the bulk viscosity 



we can calculate derivatives w.r.t. T 2 . Thus all quantities 
in (|2"5")) are uniquely determined by PHSD. We recall that 
the DQPM description of thermodynamic properties of 
1QCD results [79| as well as its implementation in PHSD 
give practically the same results |74| . 



E. Results for the bulk viscosity 

In Fig. [S] we show the bulk viscosity to entropy density 
ratio £/s as a function of temperature T of the system 
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FIG. 6: (color online) The specific sound channel (?7 + 3£/4)/s 
as a function of temperature T of the system obtained by the 
PHSD simulations in the box using the relaxation time ap- 
proximation including mean- field effects (red line+diamonds). 
It is compared with lattice QCD data from [4l|] (green circles) 
and from combining results of [1^] and [IFj] (blue squares). 



FIG. 7: (color online) The bulk to shear viscosity ratio as a 
function of temperature of the system obtained by the PHSD 
simulations in the box employing the relaxation time approx- 
imation including mean-field effects (red line+diamonds) and 
without potential effects (blue line+circles). Note the loga- 
rithmic scale in (,/eta. 



obtained by the PHSD simulations in the box employing 
the relaxation time approximation including mean-field 
(or potential) effects (red line+diamonds) and without 
potential effects (blue line+open triangles) as well as the 
available lattice QCD data from 38, 40] . Without mean- 
field effects we find an almost constant ratio (,(T)/r](T) 
(see below) which is not in line with the findings from 
the lattice. Thus the dynamical mean-fields (as incorpo- 
rated in PHSD) play a decisive role for the temperature 
dependence of the bulk viscosity ({T). 

Further quantities are of interest, in particular the spe- 
cific sound {rj + 3£/4)/s. A sound wave propagation in 
the z-direction with wavelength A = 2n/k is damped ac- 
cording to 



T 03 (t,k) cx exp 



C)kH 



2(e + p) 



(26) 



where X03 is the momentum density in the z-direction, e 
is the energy density and p is the pressure. Thus both the 
shear r\ and bulk C viscosities contribute to the damping 
of sound waves in the medium and provide a further con- 
straint on the viscosities. In Fig. [5] we present the specific 
sound channel (?y+3C /4) / s as a function of temperature T 
of the system obtained by the PHSD simulations in the 
box using the relaxation time approximation including 
mean-field effects (red line+diamonds). It is compared 
with lattice QCD results for pure SU C (3) gauge theory 
from [4l[ (green circles) and from combining results of 
37j and [38[ (blue squares). Note that the PHSD calcu- 
lations correspond to unquenched 3-flavor QCD and thus 
are not expected to match the results for the pure gauge 
theory exactly! 



Finally, in Fig. we show the bulk to shear viscos- 
ity ratio Q/r] as a function of temperature of the system 
extracted from the PHSD simulations in the box using 
the relaxation time approximation including mean-field 
(or potential) effects (red line+diamonds) and without 
potential effects (blue line+circles). Whereas an almost 
temperature independent result is obtained when dis- 
carding mean-field effects, a strong increase close to T c is 
found when including the mean-fields from PHSD. Thus 
the viscosities within PHSD give a minimum in the shear 
viscositiy T] and a strong maximum in the bulk viscosity 
C close to T c (note the logarithmic scale). 



IV. SUMMARY AND CONCLUSIONS 

We have employed the off-shell Parton-Hadron-String 
Dynamics (PHSD) approach in a finite box with periodic 
boundary conditions for the study of the shear and bulk 
viscosities as functions of temperature (or energy density) 
for dynamical infinite partonic and hadronic systems in 
equilibrium. The PHSD transport model which is based 
on a lattice QCD equation of state [7!| and accordingly 
well describes the entropy density s(T), the energy den- 
sity e(T) as well as thepressure p(T) in comparison to 
the 1QCD results [13, [H, [74|. We have employed the 
Kubo formalism as well as the relaxation time approx- 
imation to calculate the shear viscosity r/(T). We find 
that both methods provide very similar results for the 
ratio rj/ s with a minimum close to the critical tempera- 
ture T c while approaching the perturbative QCD (pQCD) 
limit at higher temperatures. For T < T c , i.e. in the 
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hadronic phase, the ratio r//s rises fast with decreas- 
ing temperature due to a lower interaction rate of the 
hadronic system and a significantly smaller number of 
degrees-of- freedom (or entropy density). Our results are, 
furthermore, also in almost quantitative agreement with 
the ratio ri(T) / s(T) from the virial expansion approach 
in Ref. [2J] as well as with 1QCD data for the pure gauge 
sector. 

We have, furthermore, evaluated the bulk viscosity 
C(T) in the relaxation time approach and focused on the 
effects of mean fields (or potentials) in the partonic phase. 
Here we find a significant rise of the ratio ((T)/s(T) in 
the vicinity of the critical temperature T c due to the 
scalar mean- fields from PHSD. The result for this ratio 
is in line with that from 1QCD calculations. Addition- 
ally, the specific sound (rj + 3£/4)/s(T) has been calcu- 
lated and presents a non-trivial temperature dependence; 
the absolute value for this combination of the shear and 
bulk viscosities is in an approximate agreement with the 
lattice gauge theory. Furthermore, the ratio ((T)/r](T) 
within the PHSD calculations shows a strong maximum 
close to T c , which has to be attributed to mean- field (or 
potential) effects that in PHSD are encoded in the in- 
frared enhancement of the 'resummed' coupling g{T). 
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